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Abstract 

Fission plays an important role in the r-process which is responsible not only for the 
yields of transuranium isotopes, but may have a strong influence on the formation 
of the majority of heavy nuclei due to fission recycling. We present calculations of 
beta-delayed and neutron-induced fission rates, taking into account different fission 
barriers predictions and mass formulae. It is shown that an increase of fission barriers 
results naturally in a reduction of fission rates, but that nevertheless fission leads to 
the termination of the r-process. Furthermore, it is discussed that the probability of 
triple fission could be high for A > 260 and have an effect on the formation of the 
abundances of heavy nuclei. Fission after beta-delayed neutron emission is discussed 
as well as difi^erent aspects of the influence of fission upon r-process calculations. 
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1 Introduction 



Fission is an important process which is responsible not only for the yields of 
the transuranium isotopes in the r-process calculation and for termination of 
the r-process, but also for the formation of the majority of heavy nuclei due 
to cycling of the r-process from transuranium region into the region of A;^130 
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(where fission products are involved again in the r- process nucleosynthesis). 
For a consistent treatment of the r-process, we have to consider fission barriers 
based on the same mass formulae that are used for the calculation of other 
reaction rates and decay properties. We have to evaluate all fission branches: 
neutron-induced fission, beta-delayed fission and spontaneous fission. In ad- 
dition, one should analyze the importance of every separate fission mode for 
the termination of the r-process, as well as nuclear abundance yields and the 
related age determination from the chronometer nuclei abundances. 

The fission barriers, used in the past for fission rate calculations [1] were prob- 
ably underestimated. Recent calculations [2,3] predict higher fission barriers, 
sometimes significantly. Therefore, a recalculation of fission rates is important 
as well as an evaluation of their sensitivity on nuclear data. We consider all 
fission channels, but first discuss beta-delayed and neutron-induced fission. 

While the study of fission in the r-process has a long history [4] , its main focus 
was related to cosmochronometers [5,6,7]. Only chronometer abundances are 
affected if the neutron freeze-out occurs when the r-process (in the envelope of 
a supernova) approaches the transuranium region and yields of fission products 
play small role for the abundances of all light r-process nuclei. For high neutron 
densities existing over long periods, fission cycling to lighter nuclei becomes 
important. For the latter case, considering only instantaneous fission is a first 
step towards a consistent treatment as performed in a number of recent studies: 
(a) Pfission = Psf = I for aU A > 240 [8]; (b) Pf.ssion = Psf = I for 
all A > 256 [6]; (c) for A = 260 [9]. In all these cases only a symmetric 
mass distribution was considered, while other calculations already included 
asymmetric mass distributions [4,7]. 

In the first extended calculation of beta-delayed fission for the majority of r- 
process nuclei [10], it turned out that the fission probabihty could reach 100% 
and beta-delayed fission affects strongly the transuranium yields, especially 
the yields of cosmochronometer nuclei. Their yields are used to determine the 
stellar age or the age of our galaxy as the lower limit on the age of the uni- 
verse. For r-process conditions with long durations of high neutron densities, 
neutron-induced fission could also play an important role. That can be seen 
from the comparison of neutron- induced fission rates for uranium isotopes [1 1] 
employing the old barrier estimates [1] and extended calculations for a variety 
of fission barrier models [12]. The availability of new extended fission barrier 
calculations [2,13], requires that all the fission rates need to be reevaluated for 
the inclusion in r-process calculations. 

The latest observations of metal-poor stars emphasize also the importance of 
the mass region 110 < A < 130 for the understanding of the role of fission in 
the r-process. One notices that these nuclei are underabundant and the solar 
system r-process must have at least two components, the one dominating for 
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A<130 and the other one dominating for A>130 [14]. The observation of the 
second component in low mctalhcity stars and the observed abundances for 
A<130 might give strong clues to fission properties, if the explanation of the 
second component is related to an r-process site with a high neutron supply 
which leads to fission cycling. This should be investigated as the site of the 
r-process is still uncertain. 

Here we want to emphasize that, besides neutron-induced fission and beta- 
delayed fission rates, the mass-distribution of fission yields is important as 
well as it was discussed earher [15]. Although the mass distribution of fission 
products is known rather well for long-lived transactinide nuclei, the mass 
distribution after fission of neutron-rich heavy elements is an open question. 
The r-process yields in the mass region of 80 < A < 120 can be affected by 
the mass distribution of fissioning nuclei, and triple or ternary [16] fission can 
play an important role in some cases. 

Beta-delayed fission rates differ very strongly due to difficulties to predict reli- 
able beta-strength functions for very neutron-rich and especially for deformed 
nuclei as well as reliable fission barriers for such nuclei. 

In a number of calculations [10,17], (contrary to [18]) the values of beta-delayed 
fission probabilities for a large number of short-lived neutron-rich heavy nu- 
clei can attain 100%, when applying the fission barriers of [1], (especially 
for transuranium nuclei) independent of beta-strength functions were used 
(see Fig. 1, Fig. 2). The beta-strength function used in [10] was based on 
Tamm-Dancoff approximation for the Gamow- Teller residual interaction. An- 
other calculation [17] used the proton- neutron Quasiparticle Random Phase 
Approximation (pn-QRPA) [19] with a shcmatic Gamow- Teller residual inter- 
action. Three sets of calculations with 3 different mass formulae [1,20,22] were 
considered taking into account an implicit treatment of nuclear deformation. 

Beta-delayed fission calculations for a number of selected nuclei were also per- 
formed [21] with a schematic consideration of deformations as in [17], to mix a 
basis of deformed nuclear states via the GT residual interaction. Fig. 1 shows 
some values of Pp^f with different parameterization of the /3-strength-function 
(squares) and artificial overestimation (triangles) and underestimation (cir- 
cles) of fission barriers and values from existing data sets [10,17] for the iso- 
topes of Pa (Z=91). Both of the existing data sets [10,17] strongly depend 
on the /3-strength function, mass models and fission barrier heights, and their 
application to r-process models gave quite different results [15]. 

The degree of difference depends on the r-process model used. In the case 
of the classical model of the r-process, the path of the r-process lies along 
nuclei with neutron separation energy Sn ~ 2MeV, and the r-process stops 
at mass numbers ARi260. Both of the data sets [10,17] yield the maximum 
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Z=91 




Fig. 1. The calculated probabilities for beta-delayed fission Pbdf for Pa- isotopes 
from different calculations using approach of Tamm-Dancoff [10] (line), QRPA [17] 
(dashed line) and RPA approximation with some artificial parametrization [21] (cir- 
cles, squares and triangles). 

values of P^j for these nuclei. That is why the previously used approximation 
of 100% instantaneous fission in the vicinity of A~260 [6], [9], [11] was a fairly 
good approximation. But when considering extreme r-process paths in the 
framework of the classical model {Sn ~ 3 MeV or S'„ ~ 1 MeV) or with a 
calculation of the heavy element abundances in a very high neutron density 
environment, the difference in the results with utilization of different beta- 
delayed fission rates can be large. 

It has to be noted that for very neutron rich nuclei the possibility of fission 
after delayed neutron emission could be rather high and, that in principle, such 
a fission mode also should be calculated. In Fig. 2, different calculations of the 
beta-delayed fission probabilities [10,17] with the same fission barriers [1] are 
compared, including the rates of fission after beta-delayed neutron emission. 
The strength-function dependence and the high values of Pnisdf are clearly 
seen, especially for heavy Cf-isotopes. 

The details of the statistical approach will be introduced in Section 2, including 
a discussion of the beta-strength function and probability of fission after beta- 
delayed neutron emission. A comparison of different fission barriers and mass 
predictions used in present calculations is shown in Section 3. In Section 4 
we present the results. The impact analysis and a comparison of the different 
decay modes in fission is discussed in subsection 4.1 with evaluations based on 
older mass predictions [22] and fission barriers [1] , often used in astrophysical 
applications. In subsection 4.2, we discuss the fission rates for more recent 
nuclear physics input from the Thomas-Fermi [2] and ETFSI [13] models and 
discuss the changes. In Sections 5 and 6, the possible consequences of fission 
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Fig. 2. Fission probabilities Pp^f (thin lines) and Ppndf (bold lines) for U and Cf 
isotopes for different strength-functions: Tamm-Dankoff [10] (lines) and pn-QRPA 
[17] (dashed lines). Mass predictions [22] and fission barriers of ref. [1] were used. 

rates are considered: cycling of the r-process due to fission into the region of 
Af»130, and the infiuence of triple fission on the yields of medium mass nuclei. 



2 The statistical approach to fission rates 



As in previous approaches [10,17], we have applied the statistical Hauser- 
Feshbach formalism for the calculation of fission rates. It has been shown [23] 
that the statistical model is very well applicable for the astrophysical neutron- 
induced rate calculations, as long as there is a high density of excited states, 
which is the case for heavy nuclei. Of course, just near the neutron drip-line the 
systematic errors of the approach in neutron-induced fission rate calculations 
can rise, underlining that reliable mass predictions are absolutely necessary for 
r-process applications far from stability. Early r-process calculations made use 
of the mass predictions by Hilf et al. [22] and the fission barriers of ref. [1] for 
fission rate evaluations. For a consistent treatment of fission rates, however, 
the neutron separation energies, reaction Q-values and fission barrier heights 
should be used from the same mass model. 



2. 1 Neutron- Induced fission and beta-delayed fission 



The cross-section for a neutron- induced reaction i^{n, /) from the target ground 
state i° with center of mass energy Ein and reduced mass is given by 

^ . 7r;.7(2//A) ^^^^ , ^, T°(^,J^£;°,J°,7^°)^/(£;,J-) 

""^^^-^ = (2jo+i).(2j„+i) i^}'' + tue;!^) • 
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The transition coefficient Tf{E, J'^) includes the sum over all possible final 
states. Since the work of Strutinski [24] fission has been generally described 
within the framework of double-humped fission barriers. Similar to previous 
work [10] we followed this approximation of a two-hump barrier. The cal- 
culation of the fission probabilities was performed in the complete damping 
approximation which averages over transmission resonances, assuming that 
levels in the second minimum are equally spaced. 

Then beta-delayed fission rate can be expressed by 



Q/3 

A/= lY.m.J.^^)-^^%^JE (2) 



where the beta-feeds Pi{E, Ji,7ri) are expressed via reduced transition proba- 
bilities B{E, JijTTi), taking into account a Gauss spreading of states: 



Ji^TTi) = -^^^exp[-{E - E,y/2a^]B{E, J^,^:^)f{Qp - E) (3) 



with integrated Fermi function 

Eq 

f{Qp -E)^l F{Z, e)e^{e'-l){Eo - efde. 
1 

Tf{E, Ji, "Kij/Tioi denotes the fission probability of the compound nucleus: 

Tf{E,Ji,7ri) _ ^ ^ Tn{E,Ji,7ri)+T^{E,Ji,7ri) y ^ 
Ttot{E,Ji,7ri) \ %s{E,Ji,Tri) J 



2 . Tn{E, Jj, TTj) + Ty{E, Jj, TTj) _ ^^^Yi[ ^a{E, Ji,ni) +Tb{E, Ji,ni) ^^^_^ 

Teff{E, Ji,7li) 2 



where 

Ta{E, Ji, TTj) ■ TeiE, Ji, Tli) 



TA{E,J„7r,)+TB{E,Ji,7ri) ' 
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with 



E 

Ta,b{E, J, tt) = j pa,b{^, J, 7r)THw{E - Ea,b - e - h%l + l)/2e, hujA,B)de . (5) 



Thw denotes the HiU- Wheeler [25] transmission coefficient through a parabohc 
barrier 

and the available energy is reduced by the rotational energy with the moment 
of inertia 9 deduced from the irrotational flow model [26]. The fission barrier 
heights Ea^b will be discussed in Section 3. Pa,b define the level densities 
above the first and second saddle points, which show an enhancement over 
the level densities at ground state deformation, due to increased deformation 
and coupling to low-lying rotational excitations. This enhancement is larger 
for the first axially asymmetric barrier than for the second mass asymmetric 
barrier. As in previous calculations, we used constant enhancement factors 4 
or 2, respectively, over the level density at ground state deformation [27,28]. 

The back-shifted Fermi-gas description of the level density was improved by 
introducing an energy dependent level density parameter a as described in 
ref.[23]. Tests show a very good agreement with shell model calculations [29] 
and justify the application at and above the neutron separation energy. 

The amount of experimentally known fission barriers used in the present cal- 
culations was increased recently by approximately a factor of 2 following the 
IAEA-report [27] and recent compilation in [13]. 

The relative probabihty of beta-delayed fission over pure beta-decay is Ppdf — 
A// A where A denotes the /3-decay rate, obtained from Eq. (5) without the 
term Tf/Ttot- 

2.2 Beta- strength function 

The beta-strength functions, defined as 
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for our calculations were derived on the basis of the QRPA calculations of 
reduced transition probabilities B{E, Ji,7rj) [31] (see formula (3)). 

In our calculations of beta-delayed fission probabilities we used the new ap- 
proach for the beta-strength function [31], based on the approach [30], consid- 
ered allowed transitions and deformed single-particle model as in [32] and takes 
into account the residual interaction. The model [31] combines calculations 
within the quasi-particle RPA for the Gamow- Teller part with an empirical 
spreading of quasi-particle strength and the gross theory for the first-forbidden 
part of /3-decay. 

In total fission barriers, penetrabilities/widths of all competing branches of 
the compound nucleus after beta-decay, and the beta strength functions are 
required to beta-delayed fission. 

Pl3df = using Eq. (2) can be written as 

Qb „ 

/ S,{E)^^dE 

PMiZ,A) = ^^ (7) 

/ Sp{E)dE 


To calculate /5-decay Q-values and neutron separation energies the experi- 
mental ground-state masses were used where available, otherwise calculated 
masses (ETFSI [33], Thomas- Fermi [2] or Howard and MoUer [1]). To reduce 
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the uncertainties in the beta-strength function calculation due to uncertainties 
in individual level energies, a spreading of individual resonances were intro- 
duced above E=2 MeV. Specifically, each "spike" in the beta-strength func- 
tion above 2 MeV was transformed into a Gaussian with width 0.5. The 
calculations were based on ground state deformations which affect the energy 
levels and wave-functions that are obtained in the single-particle model. The 
ground-state deformations were calculated in the FRDM mass model [34]. In 
some cases the first-forbidden strength calculated in the gross theory [35] was 
included. 

The example of beta-strength function for the decay of ^^'^JJ is shown in 



2.3 Cascade Model for secondary fission 

All heavy nuclei near the r- process path are very neutron rich. Therefore the 
neutron separation energies for the daughter nuclei, produced in beta-decay, 
are very low. As a result, the probability of beta-delayed neutron emission can 
be high, frequently higher than probability of beta-delayed fission, (depending 
on fission barrier value). Therefore the excited states of a daughter nucleus 
produced by the beta-decay of neutron rich transuranium nuclei, mainly decay 
by neutron emission or fission. If neutron emission is the dominant decay mode 
it can lead to excited states that lie above the respective fission barrier and 
fission could occur in the second decay-step. This probability of secondary 
fission - Ppdnf-i i-G-; fission of the compound nucleus that is produced by neutron 
emission of the excited state of the daughter nucleus, will contribute to the 
total fission probabihty P|^j and might even enhance it significantly: 



The previous calculations [10,17] showed that P^^n/- values can in certain cases 
be larger than P/^^^-values, especially for neutron-rich isotopes of the Am-Cf 
region. In Fig. 2 some of the former results were shown for the mass relations 
of Hilf et al. [22] and fission barriers of ref. [1] (Fig. 4). 

In order to evaluate the role of secondary fission on the base of new data 
this process can be investigated by iterating the SMOKER-code: if neutron 
emission leads to an excited level of the residual nucleus, we again calculate the 
branching ratios for the subsequent decay with the statistical model including 
the fission channel. Thus the probability for fission in the second decay step is 
essentially given by multiplication of the branching ratios for neutron emission 
and fission of the residual nucleus. Note that this is done by keeping track of 



(Fig. 3). 




(8) 
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Fig. 4. The neutron separation energies (line) and fission barriers (dotted line) for 
U (left) and Cf (right) isotopes on the base of Howard and MoUer (1980) approach. 

energies, angular momenta and parities of the excited states, and the final 
probability for (/3~, n/)-decays is obtained by summing up the contributions 
from all quantum mechanically allowed partial waves. 

Since the amount of excitation energy in the system is significantly reduced 
with every particle emission and eventually drops below the particle emission 
thresholds or fission barrier, we found that the cascade, could practically be 
terminated after two emission steps, contrary to neutrino-induced fission [37] 
when some subsequent emissions of different particles can occur. Thus, the 
residual daughter nuclei after two iterations were assumed to be in their ground 
state. 

The significant values of fission probability after beta-delayed neutron emission 
(second step of cascade fission) - P^dnf can be reached only for the beta-decay 
with small beta-delayed fission probability P/^df (first step of the cascade) and 
rather high beta-delayed neutron emission probability Pf^dn of mother nucleus. 
Especially high values of fission probabilities after neutron emission (second 
step of the cascade) can be achieved during decay of even-even mother nuclei 
(A,Z), when neutron separation energy in final nucleus (Z+1,A-1) significantly 
less than in daughter nucleus (Z+1,A) (See decay of U with A-266, 268, 270, 
Table 1 in Section 4.2.2). 



3 Survey of fission barrier and mass predictions 

As astrophysical applications include unstable nuclei not yet attained by ex- 
periment, all physical quantities have to be obtained from theoretical pre- 
dictions. The quantities, which need to be employed for the calculation of 
neutron-induced and beta-delayed fission include the reaction Q-values and 
the fission barrier heights. For stable nuclei, this information can come from 
the experiment, but for unstable nuclei one has to make use of a mass formula 
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and theoretical predictions of fission barrier heights. 

In addition to the Q-values and fission barriers involved, neutron emission 
and beta-decay require the knowledge of neutron optical potentials and beta 
strength functions. The accuracy of all these quantities strongly infiuences 
the probabilities of beta-delayed fission [11,18,21,36]. Several predictions exist 
so far, but their results differ strongly due to uncertainties for beta-strength 
functions and fission barriers of very neutron-rich and especially for deformed 
nuclei. 

Fig. 2 and Fig. 4 show the connection between fission barriers and neutron 
separation energies on the one hand and fission probabilities, on the other 
hand, when Howard and MoUer predictions were used. When fission barriers 
are of the order of neutron separation energy or become smaller than Sn, Ppa} 
increased and reached high values. 

We also considered two different up-to-date fission barrier predictions, devel- 
oped during recent years: one being based on the Thomas-Fermi model [2] 
and another using the Extended Thomas-Fermi model (ETFSI) [13]. Both the 

predictions, used in our calculations, include also the shell corrections. The 
analysis of fission barriers, predicted by these approaches, shows that the ear- 
lier fission barriers calculated by Howard and MoUer [1] are systematically 
underestimated . 

Fission barriers of the ETFSI model [13] made use of the (c,h)-parametrization 
[38] (see detailed explanation also in [28]). According to the ETF-predictions, 
fission barrier values increase with increasing N and that is why ETFSi pre- 
dictions [13] show significantly greater fission barrier values than the MS cal- 
culations [2] due to several reasons, mentioned in [3,13]. The shell corrections 
due to the second term of total energy based on Strutinsky Integral can even 
increase the values for nuclei with magic numbers. 

The accuracy of fission barrier predictions for large neutron excesses is hard 
to evaluate, and its effect on r-process results has to be studied. 

The second approach [2] is based on self-consistent, semiclassical mean-field 
solution of the problem of self-bound nucleons interacting by suitably adjusted 
effective forces. This modification of Thomas-Fermi model of nuclei has been 
applied to the calculation of nuclear masses and fission barriers. To predict 
fission barriers in a range of N and Z of interest a simple algebraic expres- 
sion was proposed. Predicted fission barriers include also the shell corrections 
(shell-correction method is described in [34]). The comparison with known 
experimental data shows the good agreement with calculations and one hopes 
the extrapolations of fission barriers are also reliable. 

In Fig. 5, the neutron separation energies, as well as fission barriers for U and 
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Fig. 5. The neutron separation energies (thin dotted hne) and fission barriers (sohd 
dashed hne)of daughter nuclei under beta-decay of U (left) and Cf (right) isotopes 
on the base of different barrier predictions [2] (upper row) and [13] (low row). 



Cf isotopes, calculated from considered mass relations, are shown. As it can 
be seen from the figure (the common features are the same for all actinides), 
the change of barrier heights as a function of neutron number is very similar 
for the different approaches, but differs strongly, especially for N ^ 184, where 
Mamdouh et al. [13] show the absence of deformation and increasing barriers 
in excess of 10 MeV. The older mass formula of Hilf et al. [22], often used 
for astrophysical applications, predicts significantly smaller neutron numbers 
for the position of the neutron drip line than refs [1,2]. This has important 
implications, because in self-consistent calculations the r-process path can 
pass nuclei with larger neutron excess. For the region of interest, A > 250, the 
differences in Sn can be larger than 1 MeV. 

While the agreement for experimentally known isotopes is quite good, the dis- 
crepancies of the approaches in the extrapolated region can achieve some MeV. 
These uncertainties can affect strongly the fission rates and the r-process. 
Therefore self-consistent calculations, predicting neutron separation energies, 
barriers and other nuclear properties on the same basis of the same hopefully 
reliable mass model are needed. 
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Fig. 6. The values of (n,7), (n,f) and /3-decay rates (in s ^) for U (left) and Cf 
isotopes. 

The preceding discussion (see also [48]) leads us to the conclusion that the 
fission barrier predictions from [1] underestimate whereas those from [13] over- 
estimate the barrier at least in the vicinity of N ^ 184. A main problem that 
arises when utilizing the barriers of ref. [2] lies in the fact that only the max- 
imum barrier height is given while our calculations are based on the double- 
hump representation of fission [24]. In that case we estimated the height of 
the second barrier due to barrier differences from former calculations [1]: 



Bf, = Bfm: 5/. =i?/([2])-A; where A = i?;,([l]) - i?^,([l]) (9) 



Before applying the new fission rates to nucleosynthesis, we will evaluate the 
importance of individual contributions from different fission processes, i.e. 
beta-delayed fission versus neutron-induced fission and their dependence on 
different mass excess and fission barriers predictions. 



4 Results 



The approach described above was employed for the calculations of neutron- 
induced fission rates for neutron-rich isotopes up to the neutron-drip line, as 
well as beta-delayed fission rates for all isotopes of U, Pu, Am, Cm and Cf 
with beta-decay energies Qp > 0. 



4-1 The competition of different decay channels in the transuranium region 



When the r-process approaches the actinide region, fission channels appear. 
Beta-delayed fission, proposed by Berlovich and Novikov [39], competes with 
beta-delayed neutron emission and pure beta-decay. Also the well known pro- 
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Fig. 7. The map of rates ratios Xnj/^^df for pYn = 1 and Tg = 1. 
The most abundant nuclei along the r-process path, for ^ 10^^ (crosses), 
and Hn ~ 10^^ (pluses) are marked. The position of the neutron drip-line (full 
line), nuclei with neutron separation energy Sn ~ 2MeV (dashed line) and 
nuclei with beta-decay energies Qjs ~ Sn (dotted line) are denoted as well. 
These results were obtained for the mass model of Hilf et al. [22] and fission 
barriers by Howard and MoUer [1]. 



cess of neutron induced fission has to be included. This process, discussed 
earlier [11] for uranium isotopes has not yet been included in r-process calcu- 
lations. Neutron-induced fission rates can be very high for isotopes of transura- 
nium elements and exceed beta-decay rates by orders of magnitude even near 
the neutron drip-line. This is mainly the case for very high neutron densi- 
ties Un as the product < cr„/f > Un competes with beta-delayed fission rates 
Xf3df- The recycling timescale, measuring the time until the neutron capture 
on fission products leads again to the formation of actinides can be signifi- 
cantly less than the duration time of the r-process. The different rates for U 
and Cf isotopes are compared in Fig. 6. The rates in that Figure are given for 
pYn = 1 and Tg = 1, the beta-decay rates are multiplied by 10^. It is clear 
that for pYn > 10^^ and for nuclei with A > 250 neutron-induced fission is 
more important than beta-delayed fission and radiative neutron-capture. 

For the neutron-star merger scenario [41,42] or polar jets in supernova [40] 
such conditions are maintained for a long time (~ Is). However, as already 
pointed out above, the fission barriers of Howard and MoUer [1] are very likely 
underestimated. 

The calculation of rates for neutron-induced fission and beta-delayed fission 
made use of the statistical model code SMOKER, extended by including fission 
channels [11,36] as discussed in section 2. 
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A comparison of beta-delayed and neutron-induced fission rates calculated 
with fission barriers of Howard- Moller [1] is presented in Fig. 7, where circles 
and squares of different sizes indicate different ratios of rates. One can see 
that for Z < 87 the beta-delayed fission channel is the dominating one. For 
Z = 88 — 90 neutron-induced and beta-delayed fission are in competition, 
and for Z > 91 neutron-induced fission dominates. The open squares indicate 
nuclei for which the beta-delayed fission rate is equal to zero. For the neutron 
rich nuclei {N > 160) this is due to the small neutron separation energies {Sn < 
Bf) where the decay of a daughter nucleus occurs via neutron emission rather 
than fission (because the majority of beta-decay transitions goes through low 
lying states). The most abundant isotopes with Y(A,Z)> 1% Y(Z) (crosses) 
are shown for the early stages for the r-process (when n„ > 10^^) and the time 
when the neutron density decreases to less than 10^° (pluses). 

Neutron-induced fission should be of significant importance for nuclei with the 
fission barrier -B/ ~ Sn or smaller and when the density of free neutrons is 
rather high. The rates of induced fission are large [11], [12], [43], and even for 
higher fission barriers [2,13] this branch of fission can be important. Prelimi- 
nary calculations will be discussed below. 

As already mentioned, in the majority of the r-process calculations fission was 
included in a very simple way: assuming 100% fission at given mass number 
A or charge number Z in the r-process path [7,8,9]. These simplified models 
permitted to investigate the behavior of nuclei formation due to cycling but af- 
fect the accuracy of nuclear abundances, in part, especially cosmochronometer 
yields [7,8,9]. 

The influence of a fission model in any astrophysical r-process scenario be- 
comes important when the duration time Tr of the r-process becomes larger 
than the fission cycling time Tcyde- In this case, the r-process can successfully 
reach the transuranium region and all fission rates, beta-delayed, and neutron- 
induced as well as spontaneous fission, should be taken into account. In case 
of an artificial termination of the r-process for A = Amax, the nuclei formed 
far beyond uranium undergoing a-decay will have incorrect abundances and 
will lead to wrong yields of cosmochronometer nuclei. 

In some of the runs and in previous calculations [44] we took into account 
spontaneous flssion only for nuclei with Z>92 and N>158 (according to known 
systematics) and we did not consider neutron-induced flssion for these nuclei. 
In our calculations, even for highest fission barriers [3], the neutron-induced 
fission rates are sufficiently large to terminate or at least significantly reduce 
the build-up of higher Z r-process elements. 

Simultaneous consideration of beta-delayed fission and neutron-induced fis- 
sion in our calculation, and spontaneous fission in an other approach [45] give 
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Fig. 8. Neutron-induced fission rates for U (left) and Cf (right) isotopes for different 
fission barriers as well as for different mass formulae. 

very similar results (see Section 6), at least within the calculation accuracy. 
However, neutron-induced and beta-delayed fission should be faster and con- 
sequently more important in the r-process than spontaneous fission, which 
can play, however, an important role during cooling in the final decay of the 
r-process products. 

But for a consistent inclusion of fission in the r-process, especially for the 
calculation of the final abundances, the specific scenarios of the r-process in 
very high neutron density environments should be considered. 



4-2 Mass model dependence of fission rates 



Here we present the neutron-induced fission rates and beta-delayed fission 
rates calculated for isotopes of some transuranium elements taking into ac- 
count different up-to-date fission barriers and atomic mass predictions. 



4-2.1 Neutron-induced fission rates 

In Fig. 8, the neutron-induced fission rates A^/ derived with the barriers of 
ref. [1] and mass predictions of [22] are compared with consistent calculations 
when all mass and fission barriers predictions were made on the base of ap- 
proaches [46], [47]. The results show that for uranium isotopes with atomic 
masses A<260 the rates using Thomas-Fermi predictions [2] are also rather 
high whereas the rates calculated on the base of ETFSI predictions [13] for 
atomic masses greater 250 become too small compared to beta-decay and 
neutron capture rates. The latter can not interrupt the formation of new 
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Fig. 9. Fission probabilities Ppaf for U, Am, Cm and Cf isotopes for mass models 
of Mamdouh et al. (Dashed line) and Myers and Swiatecki (line) . 

transuranium elements in r-process nucleosynthesis. As mentioned above the 
comparison was made for very high neutron environments with pYn = 1. As 
Z increases the rate values go up and all approaches give very high values of 
neutron-induced fission rates already for Cf isotopes, with the exception of 
the ETFSI approach in the vicinity of neutron shell N=184 (Fig. 8), where 
Sn « Bf. But the heights of barriers [13] for these nuclei seem to be overes- 
timated [48]. 



4-2.2 beta-delayed fission rates 

The first extended results for beta-delayed fission probabilities (Fig. 9) for the 
majority of neutron- rich isotopes of some chemical elements show that for 
fission barriers of refs. [2], [13], which are significantly higher then previously 
used [1], the values of Pf^df are still high for very neutron rich isotopes. For 
majority of them this applies in the case when mass model [2] was used and 
for part of them when the model based on ETFSI [3] was used. 

Calculations of fission rates in the framework of the statistical model (see [10] 
and Section 2 of the present work) for both self-consistent approaches [13], 
[2] show similar results when the maximum values are reached for isotopes 
with atomic mass numbers A?a 255 — 265 . But with increasing Z the region 
of isotopes with high values of P/j^/ close to 100% expands and in case of data 
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Table 1 

Beta-delayed cascade fission probabilities (in %) for some U isotopes.^-* 
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set [2] the majority of isotopes attain close to maximum values above Z—95 
and in case of [13] above Z=98. The suppression in Pf^^f values for californium 
for the ETFSI barriers is explained by very high fission barriers in the vicinity 
of nuclei with number of neutrons close to 184. 



In Table 1, the results of the cascade calculations for beta-delayed fission are 
shown. We present the results of the calculations for the nuclei with very 
small values of direct beta-delayed fission probability P/sdf- As it was expected 
the fission probability after delayed neutron emission Ppdnf can be high. For 
some decays total P^^f which is the sum of beta-delayed fission (first step of 
cascade fission) and fission after beta-delayed neutron emission (second step of 
cascade) increases by several times. That means, that the cascade mechanism 
can be important for nuclei for which the direct (one step) beta-delayed fission 
is small and beta-delayed neutron emission is high. 



In the last two columns the results of our evaluations are presented. The first 
{PpXnf) shows the fission contribution during beta-decay of daughter nucleus. 
The second one (P^^nf) shoes the same but in assumption that compound 
nucleus does not changed during all the steps of cascade. It was calculated as 
probability of delayed fission of daughter nucleus but with reduced /3-strength 
function of the mother nucleus S^iU) where U = — as in [51]. Different 
mass and barrier predictions for U [1] and Cf [13] were used. 
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Fig. 10. ETFSI and Thomas-Fermi beta-delayed fission probabilities 
5 Discussion and outlook. 



5.1 Consequences of new rates utilization for the r-process) 



The discrepancy of fission rates due to different mass predictions is high, 
especially for the region of nuclei with N ~ 184, where different mass models 
predict different dependence of fission barriers on N. New efforts are needed 
to improve the models for better agreement, mainly in the region of nuclei 
with very big neutron excess so important for the r-process. It should be 
noted that the results of Hartrcc-Fock mass calculations [48] point out that 
ETFSI-calculation [13] are overestimated, at least in the vicinity of N=184. 
That is why we can consider the fission rates of Howard and MoUer [1] were 
underestimated, the rates based on ETFSI are overestimated. This assumption 
lets us discuss what can be expected when the fission rates are applied to the 
r-process calculations. 

Looking at Fig. 10 and Table 2, it is clear that the beta-delayed fission prob- 
abilities arc quite high for the majority of isotopes of transuranium elements, 
independent of differences in mass models and strength-function calculation. 

When moderate neutron densities are considered (n„ ~ 10^^ — 10^^ for tem- 
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perature of 1.3x10^ K), the r-process passes through the region of nuclei with 
neutron separation energies Sn ~ 2 — 3MeV. Very high fission barriers, as pre- 
dicted by some calculations [13] will not avoid an influence on the r-process, 
they only move the effective region to nuclei with atomic mass numbers A 
260, where fission rates, particularly beta-delayed fission, are high (Fig. 9, 
Fig. 10). For very high neutron densities (n„ >> 10^^), when the path of the 
r-process lies closer to neutron drip line, beta-delayed fission rates for fission 
barriers by Mamdouh et al. [13] will be small. But neutron-induced fission 
rates can be significant for such conditions and additional calculations of the 
r-process are needed to clarify the question. 



In every case of long neutron exposure (with durations of the r-process >> 
Tcycie), the r-proccss seeds capture on average more than 150 neutrons, thus 
moving high abundances to the end of r-process path, where the fission stops 
the formation of heavier nuclei. Due to fission, the fission products enter the 
r-process again at smaller mass numbers (r-process cycling). The termination 
of the r-process due to the fission occurs in the range of nuclei with charge 
numbers 92 < Z < 98. Here, Tcydc - the time of r-proccss propagation from 
the average fission products to the transuranium region is defined as : 

Zi j 

For higher fission barriers a longer duration time is encountered and heavier 
nuclei are formed. 



5.2 Possible role of triple/ternary fission in the r-process 

Utilization of the latest calculations of fission barriers [2,13] should result in a 
decrease of fission rates in comparison to former barriers [1] previously used in 
r-process calculations. This in turn should lead to an increase of the maximum 
atomic numbers of fissioning nuclei formed in the r-process in excess of Z^96. 

In comparison to the Hilf et al. mass formulae [22], recent consistent ap- 
proaches shift the neutron drip line in the direction of heavier masses. For 
example, the mass prediction of Hilf et al. [22] defines the atomic mass of the 
heaviest Cf isotope to be A=279, while the other calculations discussed above 
[2,13] predict the existence of Cf isotopes heavier than A=300. The widening 
of the r-process region can result in a shift of the r-process path into the di- 
rection of larger atomic mass numbers. In case of r-process environments with 
moderate neutron densities binary fission leads to an abundance curve down- 
ing strongly below ^4 ~ 130 due to cycling of the r-process after all hghter 
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Table 2 

Beta-delayed fission probabilities for U and Cf isotopes. 
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Mass predictions were made according to ref. [22] 



nuclei have vanished due to neutron capture. For very high neutron density 
environments, when the r-process passes partly along or very close the neu- 
tron drip-line, the mass of light fission fragment can be greater than 130. Such 
a mass distribution of fission fragments could affect the good agreement of 
calculations and observation of heavy element abundances when using binary 
fission only. 

Triple or ternary fission [16,49,50] can populate nuclei with A < 130 and 
should be considered as well. It can save the agreement with observations, but 
to define the exact contribution of triple fission into the r-process the extended 
nucleosynthesis calculations should be done. 

Ternary and quaternary fission are observed experimentally [52]. Theoretical 
evaluations of the possibility of fission into n equal fragments (ternary fission 
when n — 3) followed [53]. For example, when the fissility parameter exceeds 
0.61 (that is usual for neutron rich transuranium nuclei), ternary fission ener- 
getically becomes more preferable to binary fission. 

Of course, the probability of triple fission of long lived or experimentally known 
nuclei is small because of the high second barrier of ternary fission and long 
path to the saddle point of ternary fission in comparison with binary fission. 
The experimental evaluation of the ternary fission probability in neutron in- 
duced reactions at low energy shows that this process is at least by 3 orders 
smaller than binary fission and the formation probability of a heavy third 
fragment (A> 40) decreases below 10~^ [50] . The latest extended experiment 
using the FOBOS-detector [54] shows that fission yields for three fragments 
with Z > 20 are approximately 10~^ lower than the binary ones in the nuclear 
system formed in the -"^^A^ +^^^ Th collision. In forthcoming experiments the 
ternary yields in the decay of nuclei with larger Z'^/A will be measured [54]. 
But experiments with heavy ions [49] show that the yield of ternary fission 
fragments of stable or long-lived isotopes of Th, U and Bi can be of the same 
order as in binary fission. The masses of ternary fragments were approximately 
equal. 

Up to now there are no theoretical predictions for ternary fission of unstable 
neutron-rich nuclei involved in the r-process [16]. According to the theory of 
ternary fission, fission into three fragments occurs if reach the quasi-stable 
ternary valley [16], which is difficult to investigate experimentally because 
of the high second barrier. But with increasing atomic number and proba- 
bly neutron number, the second barrier of the ternary valley decreases, espe- 
cially when the lightest third fragment A^ is significantly heavier than alpha- 
particles, and the probability of ternary fission increases. 

All these experimental facts and theoretical evaluations let us consider that 
after beta-delayed fission of extremely neutron rich nuclei (when the energy of 
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beta-decay is about 15 MeV or more and is of the order of shell corrections), 
the probability of ternary fission of nuclei with fissility x > 0.61 can be rather 
high. 

In the extreme case when ternary fission of very neutron rich nuclei is strong, 
many intermediate nuclei can be formed in the r-process and the impact of 
such a fission process upon the yields of intermediate r-elements should be 
investigated. 

At least two extreme cases can be discussed: 

1) true ternary fission into 3 approximately equal fragments (quatro fission 
occurs when x > 0.87 and it is not considered, because the fissility values for 
neutron-rich nuclei involved in r-process lie mainly in the x=0.55-0.7 region). 

2) triple fission when the average mass of the third fragment is small: A3 ~ 
10 — 20. In this case two other fragments should have smaller masses than 
in the case of simple binary fission because 10-20 mass units can be gathered 
in the third fragment. At least one heavy fragment will be included into the 
r-process before the A — 130 peak. 



6 Summciry 

The new calculations of fission rates considered two consistent approaches. The 
calculations showed that, on a whole, the new rates based on fission barriers 
of [2,13] are lower than the ones previously used [10,12,43]. However, taking 
into account a cascade model of fission after beta-delayed neutron emission, 
the total value of fission rates can increase in some cases, especially when 
beta-delayed fission of the mother nucleus is low and primary beta-delayed 
neutron emission is high. For such nuclei the addition to delayed fission after 
neutron emission can be significant and the total value of beta-delayed fission 
can increase strongly. 

Unfortunately, the agreement of different calculations (used in different mass 
and fission barrier predictions) have not converged yet, especially for very 
neutron rich nuclei in the region with N ~ 184. 

The extended calculations of fission rates taking into account Thomas-Fermi 
[2] and ETFSI [13] fission barriers have not been finished yet. A detailed com- 
parative analysis of the influence of fission modes on the results of nucleosyn- 
thesis was made in full network calculations [43,55] using the mass relations 
[22] and fission barriers [1] , used for years in many astrophysical applications. 
An evaluation of the relative contribution of neutron-induced and /3-delayed 
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fission on the formation of r-clcments under conditions of higli neutron den- 
sities (lO^o < n„ < 10-'^°) was made. Our present calculations of fission rates 
should not change significantly our previous conclusions concerning compe- 
tition of beta-delayed fission and neutron induced fission but the resulting 
changes in r- process abundances can be significant. 

The calculation of the yields Y(A) for the range of atomic masses 120 < A < 
240 depend weakly on the model of the fission approximation. For predict- 
ing the yields of nuclei heavier than 240, especially nucleo-cosmo chronometer 
yields and abundances of intermediate nuclei with atomic mass numbers 100< 
A <120, fission should be included as precisely as possible. Nuclei fighter than 
120 can be formed if the fission of many transuranium nuclei is mainly asym- 
metric. To consider all the opportunities, a schematic model for the mass 
distribution [44], including both symmetric and asymmetric [56] fission was 
applied earlier. 

In the classical model of the r-process, the r-process path proceeds under 
moderate free neutron densities and close to nuclei with a neutron separa- 
tion energy S'„ ~ 2MeV, for which in the transuranium region the values of 
beta-delayed fission are the highest ones. According to the existing fission 
probability calculations P^df [10,17], based on old fission barrier predictions 
[1], the r-process would stop in the region of transuranium nuclei at atomic 
mass numbers A ?a 250-260 and charge Z ^90-96. In that case, the approxi- 
mation of 100% instantaneous fission in the vicinity of A~260 (for simplicity 
A=260) used in [6,9,11] would be quite good. However, new fission rates can 
change the termination point of the r-process to the region of A pa 260-280 
and Z Ri94-98. 

When considering dynamical r-process models at extremely high neutron den- 
sities, in which fission occurs strongly, the path of the r-process lies close to 
neutron drip-line, and the efficiency of beta-delayed fission depends strongly 
on fission rates for the group of nuclei close to the neutron shell N?al84. In 
particular, in the scenario realized in the model of neutron star mergers [41], 
[42], r-process nucleosynthesis proceeds mainly along the neutron drip-line [57] 
due to very high neutron densities. In this case, the role of neutron-induced 
fission will be comparatively more important than beta-delayed fission. For 
realistic fission barriers [2,13] and mass predictions [46,47], nuclei with masses 
as large as A~300 could be formed. In this case the predicted abundance curve 
for nuclei below A ~ 130 should also consider triple (ternary) fission. 

The probability of triple (or ternary) fission [49,50] is highly uncertain but 
for nuclei with A > 260 near the neutron drip-line, fission into 3 fragments 
is energetically possible and could occur. Present [50] and forthcoming [58] 
experiments should clarify the probabilities of such a processes. 
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For the development of the r-process models a detailed analysis of fission prop- 
erties is required. Except for neutron-indTiced fission and /3-delayed fission, 
neutrino- induced fission was also proposed recently [59]. In this connection 
one should notice that neutrino rates [60] can also have great importance 
for specific astrophysical conditions. Precise fission rates and mass fragment 
distributions after fission have to be obtained and included in r-process calcu- 
lations. 

After extended calculations of both beta-delayed fission probabilities and neutron- 
induced fission rates for all nuclei involved in the r-process, we hope to clarify 
the question of the role of fission rates at different stages of nucleosynthesis 
for more specific astrophysical conditions. Due to our new results (see Fig. 10), 
it seems that fission will affect the r-process, but further detailed studies are 
ahead. 
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